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We have used neutron scattering techniques that probe time scales from 10~^^ s to 10~^ s to 
characterize the diffuse scattering and low-energy lattice dynamics in single crystals of the relaxor 
PbMgi/3Nb2/303 from 10 K to 900 K. Our study extends far below Tc = 213 K, where long-range fer- 
roelectric correlations have been reported under field-cooled conditions, and well above the nominal 
Burns temperature Td ~ 620 K, where optical measurements suggest the development of short-range 
polar correlations known as "polar nanoregions" (PNR). We observed two distinct types of diffuse 
scattering. The first is weak, relatively temperature independent, persists to at least 900 K, and 
forms bow-tie-shaped patterns in reciprocal space centered on (/lOO) Bragg peaks. We associate this 
primarily with chemical short-range order. The second is strong, temperature dependent, and forms 
butterfly-shaped patterns centered on (^00) Bragg peaks. This diffuse scattering has been attributed 
to the PNR because it responds to an electric field and vanishes near Td ~ 620 K when measured 
with thermal neutrons. Surprisingly, it vanishes at 420 K when measured with cold neutrons, which 
provide ~ 4 times superior energy resolution. That this onset temperature depends so strongly on 
the instrumental energy resolution indicates that the diffuse scattering has a quasielastic character 
and demands a reassessment of the Burns temperature Td- Neutron backscattering measurements 
made with 300 times better energy resolution confirm the onset temperature of 420 ± 20 K. The 
energy width of the diffuse scattering is resolution limited, indicating that the PNR are static on 
timescales of at least 2 ns between 420 K and 10 K. Transverse acoustic (TA) phonon lifetimes, which 
are known to decrease dramatically for wave vectors q « 0.2 A"^ and Tc < T < Td, are temper- 
ature independent up to 900 K for q close to the zone center. This motivates a physical picture in 
which sufficiently long-wavelength TA phonons average over the PNR; only those TA phonons hav- 
ing wavelengths comparable to the size of the PNR are affected. Finally, the PMN lattice constant 
changes by less than 0.001 A below 300 K, but expands rapidly at a rate of 2.5 x 10~^ 1/K at high 
temperature. These disparate regimes of low and high thermal expansion bracket the revised value 
of Td, which suggests the anomalous thermal expansion results from the condensation of static PNR. 

PACS numbers: 77.84.Dy, 77.65.-j, 77.80.Bh 



I. INTRODUCTION 



The concept of nanometer-scale regions of polarization, 
randomly embedded within a non-polar cubic matrix, has 
become central to attempts to explain the remarkable 
physical properties of relaxors such as PbMgi/3Nb2/3 03 
(PMN) and PhZni/gNhz/aOa {PZK)}-^-^-^ The existence 
of these so-called "polar nanoregions" (PNR) was first 
inferred from the optic index of refraction studies of 
Burns and Dacol on PMN, PZN, and other related sys- 
tems;^ and later confirmed using many different experi- 
mental techniques including x-ray and neutron diffrac- 
tiouf^iiiSii^ ^°^Pb NMRfiS and piezoresponse force mi- 
croscopy. '^'^ Early small- angle x-ray scattering and neu- 
tron pair distribution function (PDF) measurements on 
PMN by Egami et al. cast doubt on the nano-domain 
model of relaxors. However, the recent PDF analysis 
of Jeong et al., which shows the formation of polar ionic 
shifts in PMN below ~ 650 K, and which occupy only 
one third of the total sample volume at low tempera- 



tures, provides convincing support for the existence of 
PNR»i2, Neutron inelastic scattering data published by 
Naberezhnov et al. on PMN offered the first dynamical 
evidence of PNR in the form of a prominent broaden- 
ing of the transverse acoustic (TA) mode that coincides 
with the onset of strong diffuse scattering at 650Kfiiii^ 
roughly the same temperature (« 620 K) as that reported 
by Burns and Dacol in their optical study of PMN. This 
temperature, commonly known as the Burns tempera- 
ture, and denoted by Td in the original paper, is widely 
viewed as that at which static PNR first appear. Likely 
condensing from a soft TO mode, distinctive butterfiy- 
shaped and ellipsoidal diffuse scattering intensity con- 
tours centered on (/lOO) and {hhQ) Bragg peaks, respec- 
tively, are seen below T<j in both PMNi^iiiii^ and PZN.l^ 
Similar diffuse scattering patterns were subsequently ob- 
served in solid solutions of PMN and PZN with PbTiOs 
(PMN-a;PT and PZN-a;PT); however these patterns ap- 
pear only on the Ti-poor (relaxor) side of the well-known 
morphotropic phase boundary (MPB)i-6'~l The polar na- 
ture of the strong diffuse scattering, and thus its associ- 
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FIG. 1: Top - temperature dependence of the TA phonon 
energy width Tta measured with thermal neutrons by Waki- 
moto et alw^ Middle - diffuse scattering intensity contours 
measured at (110) below (300 K) and just above (650 K) the 
Burns temperature Td using cold neutrons by Hiraka et alw^ 
Bottom - diffuse scattering temperature dependence mea- 
sured with cold neutrons (open circles) and thermal neutrons 
(dashed line) by Hiraka et al^ 



ation with the formation of PNR, was unambiguously 
established by several electric field studies of PM N^^i^^ 
and PZN-8%PT^i^i25^ all of which showed dramatic 
changes in the shape and intensity of the diffuse scatter- 
ing as a function of field strength and field orientation. 
The Burns temperature Td thus represents what is ar- 
guably the most important temperature scale in relaxors 
and is several hundred degrees Kelvin higher than the 
critical temperature Tc (which for PMN is w 210 K, but 
defined only in non-zero electric field). 

Recent studies of the neutron diffuse scattering in 
PMN by Stock et al.^ and especially those by Hlinka 
et al^ and Gvasaliya et al^^^ have proven to be ex- 
tremely important because they demonstrated the util- 
ity of cold neutron spectroscopy to the study of relax- 



ors. Generally speaking cold neutrons are ill-suited to 
lattice dynamical studies because the longer wavelengths 
necessarily limit their use to the study of comparatively 
fewer (lower-Q) Brillouin zones. On the other hand cold 
neutron spectrometers provide significantly better energy 
(TiLj) and wave vector (q) resolution than do their ther- 
mal neutron counterparts. In addition, both the (110) 
and the (100) Brillouin zones of PMN are accessible us- 
ing cold neutrons with wavelengths A « 4.26 A. The com- 
bined cold and thermal neutron study of PMN by Hiraka 
et al. exploited this fact and uncovered several major new 
findings, two of which are summarized in Fig. [TJ^* The 
first finding is that the temperature at which the strong 
diffuse scattering vanishes depends on whether the mea- 
surement is made with cold or thermal neutrons, i. e. 
the value of depends on the energy resolution of the 
spectrometer. This fact, which is illustrated in the bot- 
tom panel of Fig. [1] indicates that the diffuse scattering 
in PMN, and most likely other relaxors, contains a sub- 
stantial quasielastic component. However, a consensus 
on this issue is lacking; while the finding of Hiraka et 
al. is consistent with the study of Gvasaliya et aZ.^i^ it 
contradicts that of Hlinka et al ^^^^^ A second major find- 
ing is displayed in the middle panel of Fig. [1] where inten- 
sity contours of the diffuse scattering measured with cold 
neutrons around (110) are shown to exhibit markedly dif- 
ferent reciprocal space geometries at 300 K and 650 K. 
These data unambiguously demonstrate the presence of 
two distinct types of diffuse scattering in PMN. As the 
PNR are absent at 650 K > Td, the physical origin of 
the bow-tie-shaped diffuse scattering cross section is be- 
lieved to come primarily from the underlying chemical 
short-range order. We shall refer to this as CSRO, al- 
though another commonly-used term for this is compo- 
sitionally/chemically ordered regions (COR),'^'' 

An intriguing perspective on the Burns temperature 
in PMN is provided in the top panel of Fig. [T] where the 
TA phonon energy width Tta, which is inversely pro- 
portional to the phonon lifetime, is plotted as a func- 
tion of temperature for PMN at the scattering vector 
Q = (2, 0.2, 0). These data were measured with thermal 
neutrons by Wakimoto et al^ and are consistent with 
those of Naberezhnov et al. in that the TA mode begins to 
broaden at Td ~ 620 K, the same temperature where the 
strong diffuse scattering first appears. These data also 
show that the TA broadening achieves a maximum value 
(minimum lifetime) near 420 K, which coincides with the 
value of Td measured with cold neutrons. These data 
raise the question of how to interpret the Burns temper- 
ature Td properly in that they paint a picture, markedly 
different from that currently held, in which the PNR are 
dynamic below ~ 650 K and condense into static entities 
only at a much lower temperature of 420 K. 

The goal of our study then is to determine the intrin- 
sic value of the Burns temperature Td and to clarify its 
relationship to the diffuse scattering, lattice dynamics, 
and structure in PMN. To this end we have carried out 
extensive measurements of the neutron diffuse scattering 
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cross section from 10 K to 900 K, far below Tc and well 
above the nominal value of T^j, that probe timescales from 
10^^^ s to 10^^ s. We find that a 300-fold improvement 
in energy resolution over that used by Hiraka et al.^ 
obtained using neutron backscattering techniques repro- 
duces the same onset temperature for the diffuse scatter- 
ing; hence the intrinsic Burns temperature Td for PMN 
is 420 K. At the same time an enormous change in the 
thermal expansion is observed near 420 K, which is in- 
distinguishable from zero at low temperature. Given the 
revised value of T^, this result implies the existence of 
a direct influence of the PNR on the intrinsic structural 
properties of PMN. 

We also present new data on the effects of the PNR 
on the lattice dynamics through measurements of the 
temperature and wave vector dependence of the long- 
wavelength TA phonon energy width Tta measured near 
(110) from 25 K to 900 K. We find that TA modes with 
reduced wave vectors q <ti 0.2 exhibit the same en- 
ergy width at all temperature whereas those with q w 
0.2 A""'^ exhibit a strongly temperature-dependent broad- 
ening similar to that shown in the top panel of FiglT] 
This behavior contrasts with that observed in thermal 
neutron studies of the TO mode, which exhibits a broad- 
ening for all q < 0.2 A~^. Previous neutron scattering 
work on PMN-60%PT by Stock et al.^ a material in 
which there is no strong, temperature dependent diffuse 
scattering, and thus no polar nanoregions, found no evi- 
dence of any TA phonon broadening. In this context, our 
data lend extremely strong support to the PNR model: 
the lifetimes of TA modes with wavelengths comparable 
in size to the PNR are strongly diminished by the PNR, 
whereas long-wavelength (low q) TA phonons simply av- 
erage over the PNR and are unaffected. 



II. EXPERIMENTAL DETAILS 

The neutron scattering data presented here were ob- 
tained using the BT9 thermal neutron triple-axis spec- 
trometer, the SPINS cold neutron triple-axis spectrom- 
eter, and the cold neutron High Flux Backscattering 
Spectrometer (HFBS), all of which are located at the 
NIST Center for Neutron Research (NCNR). On BT9, 
measurements of the phonons and diffuse scattering 
were made at a fixed final (thermal) neutron energy 
Ef =14.7 meV (A = 2.36 A) using the (002) Bragg reflec- 
tion of highly-oriented pyrolytic graphite (HOPG) crys- 
tals to monochromate and analyze the incident and scat- 
tered neutron beams, respectively. Horizontal beam col- 
limations were 40'-47'-S-40'-80' (S = sample). A spe- 
cial, non-standard high q-resolution configuration was 
employed to measure the thermal expansion in which 
the (004) Bragg reflection from a perfect Ge crystal 
was used as analyzer and horizontal beam coUimations 
were tightened to 15'-47'-S-20'-40'. The choice of Ge 
was motivated by the close matching between the PMN 
(022) (1.431 A) and Ge (004) (1.414 A) d-spacings, which 



provides a significant improvement in the instrumental 
g-resolutionj^ On SPINS, which sits on the cold neu- 
tron guide NG5, the phonon and diffuse scattering mea- 
surements were made at a fixed final neutron energy 
Ef = 4.5 meV (A = 4.264 A) also using the (002) Bragg 
reflection of HOPG crystals as monochromator and ana- 
lyzer. A liquid- nitrogen cooled Be filter was located after 
the sample to remove higher order neutron wavelengths 
from the scattered beam, and horizontal beam coUima- 
tions were set to guide-80'-S-80'-80'. The resultant elas- 
tic {hu! = 0) energy resolution for the SPINS measure- 
ments was SE = 0.12 meV half-width at half-maximum 
(HWHM). 

The High-Flux Backscattering Spectrometer was used 
to look for dynamics that might be associated with 
the strong diffuse scattering below Td- This instru- 
ment uses a mechanically-driven Si(lll) monochromator 
to Doppler shift the energies of incident neutrons over 
a narrow range centered about 2.08 meV. Neutrons are 
backscattered from the monochromator and proceed to- 
wards the sample where they are scattered into a 12 m^ 
array of Si (111) crystals that serve as analyzer. These 
neutrons are then backscattered a second time by the 
analyzer, which selects the final neutron energy Ef = 
2.08 meV, into a series of detectors positioned about the 
sample. The effective angular acceptance of each detector 
is ~ 15°. The HFBS instrument is described in further 
detail elsewhere<^ The elastic energy resolution for the 
HFBS measurements described here was SE = 0.4 /leV 
(HWHM). 

Two high-quality single crystals of PMN, labeled 
PMN #4 and PMN #5, were used in this study; both 
were grown using a top-seeded solution growth tech- 
nique The crystal growth conditions were determined 
from the pseudo-binary phase diagram established for 
PMN and PbO. The PMN #4 and #5 crystals weigh 
2.7 g (0.33 cm^) and 4.8 g (0.59 cm^), respectively. At 
300 K the mosaic of each crystal measured at (220) is less 
than 0.04° full- width at half-maximum (FWHM). Loss of 
PbO, the formation of a pyrochlore phase, and the reduc- 
tion of Nb^+ are known to occur in PMN single crystals 
when subjected to high temperatures under vacuum for 
extended periods of time. This process results in a dra- 
matic blackening of the crystal, which is normally of a 
transparent gold/amber color. While dielectric measure- 
ments on such darkened crystals reportedly show little 
difference from those on unhealed samples)^ our mea- 
surements reveal a diminishment of the diffuse scattering 
intensity after sustained and repeated heating. Therefore 
experiments on the larger PMN crystal #5 were limited 
to 600 K or less, while PMN crystal #4 was used to ob- 
tain data above 600 K. 

Both samples were mounted with an [001] axis oriented 
vertically, giving access to reflections of the form {hkO). 
For the high temperature experiments, PMN crystal #4 
was wrapped in quartz wool, mounted in a niobium 
holder secured by tungsten wire, and then loaded into 
a water-cooled furnace capable of reaching temperatures 
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FIG. 2: Diffuse scattering intensity contours in PMN mea- 
sured on SPINS near (100) are sliown using a linear gray scale 
at (a) 300 K, (b) 650 K, and (c) 900 K. Diffuse scattering in- 
tensity contours measured on BT9 at 900 K near (300) using 
thermal neutrons are shown in (d). The data shown in panels 
(a) and (b) are from Hiraka et al^ Schematic diagrams of 
the diffuse scattering intensity contours below (Low-T) and 
at Td (High-T) are shown in the upper right. 



from 300 K to 1800 K. PMN crystal #5 was mounted onto 
an aluminum sample holder assembly placed inside an 
aluminum sample can, and then loaded inside the vac- 
uum space of a closed-cycle *He refrigerator that provides 
access to temperatures from 10 K to TOOK. Each sample 
has a cubic lattice spacing of a = 4.05 A at 300 K; thus 
Irlu (reciprocal lattice unit) equals 27r/a — l.SSA^^. 



III. ORIGINS OF THE DIFFUSE SCATTERING: 
PNR VERSUS CSRO 

The diffuse scattering in PMN crystal #5 was studied 
by Hiraka et alM- at and below the nominal Burns tem- 
perature Td with good energy resolution [dE k, 0.12meV 
HWHM) using cold neutrons. Data were measured in 
the (100) and (110) Brillouin zones and are represented 
schematically in the upper right portion of Fig. [51 The 
horizontal and vertical axes in this figure correspond 
to the components of the scattering vector Q = (hkO), 
which are measured in reciprocal lattice units (rlu). The 
"Low-T" (T < Td) regime is dominated by contributions 
from the PNR, where the diffuse scattering intensity con- 



tours near (100) resemble a butterfly and those near (110) 
resemble an ellipse for which the long axis is oriented 
perpendicular to Q. This is shown in Fig. [11(a) where 
diffuse scattering data near (100) at 300 K reveal an in- 
tense butterfly-shaped pattern. We note here that we 
see no evidence of any "transverse component" to the 
diffuse scattering in any of our PMN crystals like that 
reported previously by Vakhrushev et al^i- We assume 
that this component was due to the crystal imperfections 
mentioned by these authors in their PMN sample, which 
gave rise to a powder ring in their data at (100). 

The same butterfly/ellipsoidal diffuse scattering geom- 
etry was shown to persist in single crystals of PMN- 
xPT and PZN-xPT in studies by, respectively, Matsuura 
et al. and Xu et al. for compositions spanning the Ti- 
poor (relaxor) side of the morphotropic phase boundary 
(MPB).32.!^ These results also completely refute those 
of La-Orauttapong et al. who reported that the orienta- 
tion of the strong diffuse scattering varies with Ti con- 
tent in PZN-xPT and concluded that the PNR orienta- 
tion changes with doping4^i^ For Ti-rich (tetragonal) 
PMN-xPT compositions just beyond the MPB, Mat- 
suura et al. found that the strong, temperature depen- 
dent diffuse scattering vanishes and is replaced by crit- 
ical scattering;22 Matsuura et al. also found that the q- 
integrated diffuse scattering intensity increases with Ti 
content on the Ti-poor side of the MPB, peaks near 
the MPB, then drops dramatically on crossing the MPB. 
This finding is significant because it suggests that an in- 
triguing and direct correlation exists between the PNR 
and the piezoelectric coefficient CZ33, which exhibits the 
same dependence on Ti content 3 A model based on 
pancake-shaped, ferroelectric domains has been used suc- 
cessfully to fit the three-dimensional diffuse scattering 
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FIG. 3: Contour plot of the difference in diffuse scattering 
intensity near (110) measured at 550 K and 400 K. The result 
is the characteristic figure "8" pattern, observed previously 
by Vakhrushev et ali^ 
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distributions measured in PZN-xPT with high-energy x- 
raysi^ii^i^ A similar type of real-space structure has 
been proposed to explain the diffuse scattering in the 
relaxor KLia;Tai_x03.''^ On the other hand, alternative 
models explaining the same diffuse scattering distribu- 
tions have also been proposed.— 

In the "High-T" regime (T > Td) there are no PNR, 
and the associated butterfly-shaped diffuse scattering is 
no longer present. The weak diffuse scattering that re- 
mains is thus argued to originate primarily from the 
underlying chemical short range order (CSRO), which 
reflects weak correlations between the Mg^+ and Nb^+ 
cations on the i3-site of the perovskite ABOz structure. 
In this regime the shapes of the diffuse scattering con- 
tours are radically different, resembling a bow-tie in both 
Q = (/lOO) and Q = (hhO) Brillouin zones in which the 
diffuse scattering extends mainly parallel to Q. The only 
difference between the contours near (100) and (110) ap- 
pears to be in the orientation of the triangular regions 
of diffuse scattering, which point in towards (100), but 
away from (110). Data taken near (100) at 650K are 
displayed in Fig. [2](b).^* At this temperature the diffuse 
scattering intensities, shown using a linear gray scale, 
are much weaker than those of the butterfly pattern at 
300 K. We further note that the intensities increase (be- 
come darker) from left to right in panels (b) and (c), 
which corresponds to increasing Q. Given the {Q ■ u)^ 
dependence of the neutron diffuse scattering cross sec- 
tion, this intensity signature implies the presence of some 
short-range, correlated displacements u since otherwise, 
if w = 0, there would be no Q-dependence. Thus the 
weak diffuse scattering is not solely due to CSRO. 

In Fig. [3] we plot the difference between the diffuse scat- 
tering intensities measured near Q — (110) at 400 K and 
550 K. The resulting contours approximately reproduce 
the "flgure 8" pattern observed previously by Vakhru- 
shev et al.r^, and indicates that the ellipsoidal diffuse 
scattering has a strong temperature dependence. By con- 
trast, the bow-tie-shaped diffuse scattering is effectively 
subtracted out in this analysis, which conflrms that it 
has little to no temperature dependence. Thus the high- 
temperature diffuse scattering is not associated with the 
Burns temperature Td or with the formation of long- 
range polar correlations below Tc. We note that such 
a strictly incoherent treatment of the high-temperature 
(CSRO) and low-temperature (PNR) components of the 
total diffuse scattering, as described here, ignores the in- 
evitable cross terms that must exist between them. On 
the other hand, if the high-temperature scattering does 
arise primarily from CSRO, then it should be largely in- 
dependent of the ionic displacements that give rise to the 
butterfly-shaped diffuse scattering (PNR) below Td- The 
relative weakness of the high-temperature diffuse scatter- 
ing compared to that at low temperatures also suggests 
such cross terms should be weak, and this appears to be 
supported by the simple subtraction analysis presented 
in Fig. [3] in that one effectively recovers the ellipsoidal 
(not bow-tie) intensity contours. For this reason we be- 
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FIG. 4: Diffuse scattering intensity measured on BT9 at 300 K 
as a function of reduced momentum transfer relative to (100) 
(open circles) and (200) (solid circles). The dotted lines iso- 
late the diffuse scattering component from the Bragg peak. 
The inset shows the orientation of the scan relative to the 
familiar butterfly pattern measured near (200) at 100 K. 

lieve it is a reasonable flrst approximation to treat the 
two diffuse scattering components as being nearly inde- 
pendent. 

We extended the diffuse scattering measurements to 
temperatures well above Td using PMN crystal #4, which 
we reserved for very high temperature experiments. Pre- 
vious data taken on PMN and PZN crystals heated to 
1000 K revealed signiflcant evidence of sample decompo- 
sition, so we limited our measurements to 900 K. Diffuse 
scattering intensity contours at 900 K for the (100) and 
the (300) Brillouin zones are presented in Fig. [2](c) and 
[2](d), respectively. Although the diffuse scattering near 
(100) is quite weak, it is consistent with the bow-tie ge- 
ometry observed in Fig. [21(b) at 650 K. We exploited the 
dependence of the neutron diffuse scattering cross sec- 
tion to obtain higher intensity by using thermal neutrons 
to access (300). The (300) diffuse scattering intensity 
contours are shown in Fig. [l](d), where the bow-tie pat- 
tern observed at 650 K is still present at 900 K. The con- 
tours arc truncated for h > 3.02 rlu because of a mechan- 
ical limit on the maximum scattering angle available on 
BT9, but the triangular region on the low-Q side of (300) 
is clearly evident. That the bow-tie-shaped diffuse scat- 
tering persists to such high temperature provides what 
is perhaps the most convincing evidence that it arises 
mainly from CSRO. 
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FIG. 5: Diffuse scattering contours measured on BT9 near 
(200) are shown at (a) 100 K, (b) 300 K, and (c) 600 K. 



During the course of our measurements we noticed 
that the diffuse scattering in PMN crystal #4 had di- 
minished and was noticeably weaker than that in PMN 
crystal #5, which we had never exposed to temperatures 
above 650 K. We also observed a broad ring of scatter- 
ing passing directly through Q = (300), which is shown 
in Fig. [21(d). This feature was never observed prior to 
heating this crystal to 900 K. Given the length of time 
spent at high temperatures, this feature most likely cor- 
responds to a powder ring arising from partially decom- 
posed regions of PMN crystal #4, which had turned en- 
tirely black after exposure to high temperatures. These 
regions do not affect any other data presented in this pa- 
per because the 900 K measurements on PMN crystal #4 
were the last ones performed on this sample. Therefore 



the powder ring appears only in Fig.[2](d). 

It is instructive to compare these results with those 
on PMN-60%PT, a composition that lies well beyond 
the morphotropic phase boundary and undergoes a first- 
order ferroelectric transition from a cubic to a tetragonal 
phase near 540 K. An extensive study of this material 
using neutron and high-energy x-ray scattering meth- 
ods found no sign of the strong, butterfly-shaped dif- 
fuse scattering at low temperaturesi^ This result lends 
further support to our association of the strong, temper- 
ature dependent diffuse scattering with the PNR, which 
are absent in PMN-60%PT. Neutron measurements on 
PMN-60%PT do, however, reveal the presence of bow-tie 
shaped diffuse scattering intensity contours at all tem- 
peratures studied, which supports the identification of 
such diffuse scattering with chemical short-range order 
between cations on the B site of the PMN perovskite 
ABO'i structure. This picture is supported by theoreti- 
cal work^^ as well as ^^Nb NMRjA^ electron microscopyr^ 
and polarized Raman scattering^-'^ measurements. All of 
these studies suggest that there is no temperature de- 
pendence to the bow-tie shaped diffuse scattering below 
K, 1000 K, which is consistent with our results on PMN 
over the extended temperature range. 

During our study of PMN we discovered that the dif- 
fuse scattering near (200) is not as weak as previously 
believed jiii^i^'^ To confirm this finding, we made de- 
tailed measurements of the diffuse scattering intensity 
near (200) at 300 K and 100 K along a trajectory in re- 
ciprocal space that follows one wing of the butterfly in- 
tensity contour; this is shown by the dashed line in the 
inset to Fig. ID The results of the 300 K scan are com- 
pared to an identical scan measured in the (100) zone, 
both of which are shown in Fig. [D These data demon- 
strate that the (100) diffuse scattering cross section, rep- 
resented by the dotted lines passing through the open 
circles, is substantially larger than that at (200), des- 
ignated by the solid circles. This result supports the 
model of Hirota et al. in which the (unexpectedly) weak 
(200) diffuse scattering cross section observed in PMN 
and other relaxors can be explained by the presence of a 
uniform shift or displacement of the PNR relative to the 
non-polar cubic matrix along the direction of the local 
PNR polarization.^^ Indirect evidence for the existence 
of this shift has been obtained from neutron scattering 
measurements of the anisotropic response of the diffuse 
scattering in PZN-8%PT to an electric field applied along 
the [001] direction,^ 

Fig- m shows diffuse scattering intensity contours mea- 
sured on BT9 at 100 K, 300 K, and 600 K near(200); these 
data illustrate that the (200) diffuse scattering intensity 
follows the same temperature dependence as that mea- 
sured in other Brillouin zones, where the diffuse scatter- 
ing is much stronger. As the temperature is raised the 
diffuse scattering intensity decreases in the same manner 
as that previously measured and observed in the (100), 
(110), and (300) Brillouin zones. This proves that the dif- 
fuse scattering measured at (200) has the same origin as 
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FIG. 6: DifTuse scattering intensity contours in PMN mea- 
sured on the NCNR Disk Chopper Spectrometer (DCS) near 
(100) are plotted on a logarithmic gray scale at 300 K. Open 
circles represent the reciprocal space locations of the detec- 
tors of the High Flux Backscattering Spectrometer relative to 
the butterfly-shaped intensity contours. Each open circle is 
a separate detector, and the detector numbers are indicated. 
These DCS data are taken from Xu et aZ.— 



that in other zones, i. e. that it is associated with the for- 
mation of PNR. At 600 K in panel (c) one can already see 
the emergence of the bow-tie-shaped diffuse scattering 
that is otherwise obscured by the stronger PNR-related 
diffuse scattering at lower temperatures. These data are 
important because they support the mode-coupling anal- 
ysis of Stock et al^, which assumes that the diffuse scat- 
tering in PMN in the (200) and (220) Brillouin zones is 
much weaker than that in the (110) zone. Thus we em- 
phasize that while the neutron diffuse scattering cross- 
section near (200) is not zero, it is small and consistent 
with previous structure factor calculations. 



IV. DIFFUSE SCATTERING DYNAMICS: A 
REASSESSMENT OF THE BURNS 
TEMPERATURE 

The reciprocal space geometry of the strong diffuse 
scattering in PMN was first characterized using x-ray 
diffraction and is consistent with the neutron scatter- 
ing data we have presented herci ^^'^^ The energy reso- 
lution provided by x-ray diffraction {5E « 1000 meV) is 
typically much broader than that of thermal neutrons 
{6E « 1 meV) ; thus it was assumed that the strong dif- 
fuse x-ray scattering originated from low-energy, soft. 
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FIG. 7: (Color online) (a) Contour plot of the diffuse scatter- 
ing intensity measured on the HFBS as a function of energy 
transfer and temperature. Contours are shown on a linear 
intensity scale; dashed lines indicate the full- width at half- 
maximum (FWHM) of the peak linewidth at each tempera- 
ture. Data were summed over detectors 10-16 as illustrated in 
Fig.[6l Panels (b), (c), and (d) show inelastic scans at 700 K, 
300 K, and 100 K. The horizontal bar in panel (c) represents 
the instrumental FWHM elastic energy resolution {25E). 



transverse optic (TO) phonons that were captured by 
the large energy resolution] ^^'^'^ However the cold neu- 
tron data of Hiraka al. provide a much narrower elas- 
tic energy resolution of « 0.12meV HWHM and show, 
unambiguously, that the diffuse scattering cross section 
contains a component that is static on timescales of at 
least ~ 6ps below 420 K as illustrated in Fig. [TJ This 
result was subsequently confirmed on a separate PMN 
crystal by the neutron study of Gvasaliya et al. which 
employed comparable energy resolution)^ Hence the ob- 
served strong diffuse scattering cannot simply be the re- 
sult of a soft, low- lying TO phonon. The TO mode must 
condense and/or broaden sufficiently to produce the elas- 
tic diffuse scattering cross section observed by Hiraka et 
al. Such a scenario is in fact suggested by the correspond- 
ing thermal neutron data taken on BT9 using a some- 
what coarser energy resolution of w 0.50 meV HWHM. 
As shown in Fig. [1] an apparent elastic diffuse scattering 
cross section is observed up to temperatures as high as 
650 K. 

Motivated by these results, we looked for evidence of 
a dynamic component of the strong diffuse scattering 
using the NCNR High-Flux Backscattering Spectrome- 



8 



PMN, Diffuse Intensity 




200 400 600 

Temperature (Kelvin) 



FIG. 8: (a) Temperature dependence of the diffuse scattering 
measured on the HFBS; data are integrated in energy from ± 
5 neV and in Q over detectors 10-16 as illustrated in Fig. [6] 
(b) Temperature dependence of the diffuse scattering intensity 
measured on SPINS at (1.05,0.95,0) by Hiraka et al?^ The 
SPINS energy resolution of 5E^12Q ^JieV HWHM provided 
the energy integration. 



ter (HFBS), which provides an elastic energy resolution 
5E = 0.4 //eV HWHM. We oriented PMN crystal #4 in 
the (hkQ) scattering plane, which is the same geometry 
used for the triple-axis studies discussed in the previous 
sections. Fig. [5] displays the resulting locations of each of 
the 16 HFBS detectors relative to the butterfly-shaped 
diffuse scattering pattern at (100) measured previously 
using the NCNR Disk Chopper Spectrometer (DCS) 3 
From this figure it can be seen that detectors 10 through 
16 sample different parts of the wings of the strong dif- 
fuse scattering. In particular, detectors 13 and 14 lie 
close to the (100) Bragg peak. Because the instrumental 
Q-rcsolution of the HFBS is relatively poor, we checked 
for the presence of Bragg contamination in the integra- 
tion analysis by removing contributions from detectors 
13 and 14. This did not change any of the results. More- 
over the PMN study by Wakimoto et al. showed the (202) 
Bragg peak intensity changes by less than 10% between 
50 K and 300 K.^^ Therefore, in the following analysis we 



have integrated the intensity from detectors 10 through 
16. 

The energy dependence of the diffuse scattering as a 
function of temperature is illustrated in Fig. [T] Panel 
(a) shows a color contour plot of the peak lincwidth in 
energy as a function of temperature after subtraction of 
a high temperature background. This background was 
taken to be the average of the intensity measured at 
600 K and TOOK because no strong diffuse scattering is 
observed with either cold neutron triple-axis or neutron 
backscattering methods at these temperatures. The dis- 
tance between the dashed lines represents the measured 
energy width (full-width at half-maximum) as a function 
of temperature. Typical energy scans are displayed in 
panels (b), (c), and (d) at TOOK, 300K, and lOOK, re- 
spectively. We see that the peak centered a.t hui — is 
equal to the instrumental energy resolution SE — 0.4 /ieV 
HWHM at all temperatures. Based on this analysis, we 
conclude that the strong diffuse scattering is elastic on 
timescales of at least r « h/SE ^ 2 ns. Our results are 
thus consistent with the neutron-spin echo study on PMN 
of Vakhrushev et al.*^ 

Fig. [H](a) shows the temperature dependence of the 
diffuse scattering intensity integrated over ±5 ^eV. These 
data are compared to those measured on PMN crystal #5 
using the cold neutron spectrometer SPINS, which are 
plotted in panel (b). That the temperature dependences 
from the backscattering and SPINS measurements agree 
within error in spite of a 300-fold (120 /ieV/0.4 /icV) dif- 
ference in energy resolution proves that static, short- 
range polar order first appears at much lower temper- 
atures than has been understood from previous data 
measured with much coarser energy resolution. These 
data thus demand a revised value for the Burns temper- 
ature: Td = 420±20K. We mention, however, that if the 
diffuse scattering energy width obeys an Arrhenius or 
some power law form, which has been suggested for spin 
glasses in Ref. [H, then a more detailed analysis based on 
data taken closer to the onset of the diffuse scattering at 
420 K will be required to confirm (or reject) the existence 
of such alternative dynamic contributions to the diffuse 
scattering energy width. 

To gain a better understanding of the apparent 
quasielastic nature of the diffuse scattering, we examined 
the temperature dependence of the low-energy transverse 
acoustic (TA) modes in greater detail. In particular 
we focused on measurements of the temperature depen- 
dent broadening of the transverse acoustic phonon over a 
range of reduced wave vectors q that approach the zone 
center [q — 0). These data are discussed in the following 
section. 



V. TA AND TO PHONONS: EFFECTS OF THE 

PNR 

An extensive series of inelastic measurements were 
made on both PMN crystals #4 and #5 in the (110) 
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Brillouin zone using the SPINS spectrometer in an effort 
to map out the temperature dependence of the TA mode 
for reduced wave vectors q ^ q^^f and q w g^,/, where 
Qwf = 0.14 rlu ~ 0.2 A^^ is the wave vector associated 
with the so-called "waterfall" effect, and below which 
the long-wavelength, soft TO modes in PMN, PZN, and 
PZN-8%PT are observed to broaden markedly at tem- 
peratures below Td.— '^'^'^1^ This zone was chosen be- 
cause the TA phonon dynamical structure factor for (110) 
is much larger than that for (100). The TA phonon en- 
ergy lineshapes were studied at two very different values 
of the reduced wave vector q, measured relative to the 
(110) zone center, for temperatures between 100 K and 
900 K. These data are presented in Fig. [51 The data 
shown on the left-hand side of this figure correspond to 
Q = (1.035,0.965,0) or g = ^/2 • 0.35 rlu = O.OSrlu = 
0.077 A^^, which is less than half of q^^f. These data 
show that the TA lineshape at this small wave vector 
remains sharp and well-defined at all temperatures and 
has an intrinsic energy width that is larger than the in- 
strumental resolution (shown by the small, solid horizon- 
tal bar at 500 K). By contrast, the data on the right- 
hand side of Fig (HI correspond to Q = (1.1,0.9,0) or 
q = V2- O.lOrlu = 0.14rlu = 0.22 which is nearly 
equal to qwf - In this case it is quite evident that the TA 
lineshape broadens dramatically, becoming increasingly 
ill-defined below Td, especially at 300 K, but then sharp- 
ens at lower temperature, e. g. at 100 K. This behavior 
is the same as that observed for the soft, zone-center TO 
mode in PMN by Wakimoto et al. measured at (200)31 
In addition to these cold neutron measurements of the 
TA mode using SPINS, other measurements were made 
using thermal neutrons on BT9 to characterize the TO 
mode lineshapes for reduced wave vectors q « q^jf from 
300 K to 900 K. In these experiments data were taken us- 
ing a water-cooled furnace with PMN crystal #4 in both 
the (300) and (100) Brillouin zones for q — q^f — 0.14 rlu 
and are shown in FiglTOl As can be seen on the left- 
hand side of this figure, both the TA and TO modes 
are well-defined at 900 K. The scattering intensity below 
1-2 meV increases sharply because of the comparatively 
coarser energy resolution intrinsic to BT9 (~ 0.5 meV 
HWHM), which uses thermal neutrons, compared to that 
on SPINS. Even so, the TA mode is easily seen at both 
900 K and 700 K. However at 500 K the TA mode is less 
well-defined. This occurs in part because the mode has 
broadened, but also because the low-energy scattering 
has increased, which results from the onset of the strong 
diffuse scattering due to the PNR. At 300 K the TA mode 
is nearly indistinguishable from the sloping background, 
and the TO mode has become significantly damped. One 
also sees a gradual softening of the TO mode from 900 K 
down to 300 K that is of order 1-2 meV. The same data 
were taken in the (100) zone, which appear in the right- 
hand portion of Fig. [TOl Essentially the same trends are 
observed in this zone, with the TA mode again becoming 
almost indistinguishable from background at 300 K. One 
difference is that the TA mode is better separated at all 
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FIG. 9: TA phonon lineshapes in PMN measured on SPINS 
with cold neutrons at Q = (1.035,0.965,0) and (1.1,0.9,0) 
from 100 K to 900 K. 



temperatures from the steep increase in scattering seen 
at low energies. These data are consistent with the con- 
jecture of Stock et al^ that the low-energy TA phonon 
is coupled to the diffuse scattering centered around the 
elastic position. The inset on the top panel on the right- 
hand side displays the full TA-TO phonon spectrum at 
900 K out to 20meV. 

The square of the energies of both the TA and TO 
modes presented in the two previous figures are plotted 
versus temperature in Fig. 1111 This is done to draw at- 
tention to the similarity between the low-energy lattice 
dynamics of PMN and the behavior expected for a con- 
ventional ferroelectric, for which the soft mode frequency 
squared varies linearly with (T — Tc). Data for the 
zone-center soft mode measured at (200) and (220) have 
been taken from Wakimoto et al. and Stock et al. and 
added to the top of Fig. [Tl] for ease of comparison i2L2S 
The data taken by Wakimoto et al. are based on en- 
ergy scans measured at constant-Q at the zone center, 
whereas Stock et al. determined the zone-center TO en- 
ergies by extrapolating from data obtained at non-zero q. 
The extrapolation technique permits phonon energies to 
be extracted in the temperature range where the zone- 
center TO mode is heavily damped. These values for the 
zone-center TO phonon energy have been confirmed by 
infrared spectroscopy.®^ What this figure immediately re- 
veals, then, is a surprising ferroelectric character of the 
zone center TO mode above Td and below Tc, and a corre- 
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FIG. 10: TA and TO phonon lineshapes in PMN measured 
on BT9 with thermal neutrons near the waterfaU wave vector 
at Q = (3, -0.14,0) and (1,-0.14,0) from 300 K to 900 K. 



spending change in the behavior of the TA mode energy 
that is bracketed by the same two temperature scales, 
which are denoted by the two vertical dash-dotted Hnes. 
We note that no such change in the acoustic phonons 
was observed in PMN-60%PT where PNR and the asso- 
ciated strong diffuse scattering are absent; therefore the 
softening of the TA mode is directly related to the de- 
velopment of PNR.~ A further interesting feature is the 
minimum that is present in all of these data near 420 K, 
which is the same temperature at which the strong diffuse 
scattering first appears when measured with high energy 
resolution, namely the revised value of the Burns tem- 
perature. Therefore, the onset of the diffuse scattering is 
directly associated with the softening of the TO mode; 
this is yet further evidence that it is associated with the 
formation of static, short-range polar correlations. 

It is very important to note that the square of the 
TA phonon energy measured at Q=(l. 1,0.9,0), which 
corresponds to q 0.14 rlu, shows a much more pro- 
nounced minimum at 420 K than does that measured at 
(5=(1. 035, 0.965,0), which corresponds to g ~ 0.05 rlu. 
This shows that long-wavelength TA phonons exhibit a 
much weaker response to the formation of static short- 
range, polar correlations. This can be understood in 
terms of a simple physical picture in which those phonons 
with wavelengths comparable to the size of the PNR are 
strongly affected (damped) by the PNR whereas longer 
wavelength phonons simply average over the PNR and 
are thus not affected by the presence of static, short- 
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FIG. 11: The square of the phonon energy is plotted ver- 
sus temperature for the zone-center TO modes measured at 
(200) (open diamonds) and (220) (solid diamonds), taken 
from Wakimoto et al. and Stock et alM^, as well as for 
three TA modes measured at (1.035,0.965,0) (open triangles), 
(1.1,0.9,0) (open circles), and (3,-0.14,0) (solid circles) . Linear 
behavior is consistent with that expected for a conventional 
ferroelectric soft mode. 



range polar correlations. This idea was previously pro- 
posed to explain the anomalous damping of the TO mode 
for wave vectors q < q^f near the zone center.^- However, 
no strong diffuse scattering is seen in PMN-60%PT and 
thus no PNR are present, even though the anomalous 
TO broadening is still observed; hence this TO broad- 
ening, which gives rise to the waterfall effect, cannot be 
associated with the presence of static, short-range po- 
lar correlations. On the other hand, the idea that the 
acoustic phonons are affected by PNR is confirmed by 
the absence of any such acoustic phonon broadening in 
PMN-60%PT. Thus PNR have a significant effect on the 
low-energy acoustic phonons over a limited range of re- 
duced wave vectors that may be related to the size of 
the PNR. In light of the diffuse and inelastic scattering 
data that have been analyzed so far, we now turn to the 
detailed measurement of the thermal expansion in PMN. 
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FIG. 12: Lattice strain of single crystal PMN derived from the 
(220) Bragg reflection measured from 30 K to 580 K on BT9. 
Data were obtained on heating using a perfect Ge(004) ana- 
lyzer, an incident neutron energy Ei = 14.7 meV, and horizon- 
tal beam coUimations of 15'-47'-S-20'-40' to obtain extremely 
high q-resolution. X-ray data from Dkhil et al. are indicated 
by the open squares for comparison.—. 



VI. THERMAL EXPANSION: INVAR-LIKE 
EFFECT BELOW Ta 

A variety of experimental techniques have been used to 
measure the thermal expansion in both single crystal and 
powder samples of PMN, which includes x-ray diffraction, 
laser dilatometry, and neutron diffraction.—:^— Here 
we present extremely high g-resolution neutron measure- 
ments (Sq w 0.0005 rlu = 0.0008 A~i HWHM) of the cu- 
bic unit cell lattice parameter a of PMN crystal #5 using 
the special triple-axis configuration described in the ex- 
perimental section. This configuration, which employs a 
perfect single crystal of Ge as analyzer, also provides an 
exceptional energy resolution 6E « 5 — 10/ieV HWHM. 
The resulting data are plotted in Fig. [12] in terms of the 
strain a/aQ—l, where ag is the lattice parameter at 200 K. 
Measurements were made on heating from 30 K to 580 K 
using a closed-cycle *He refrigerator after having first 
cooled to 30 K. X-ray data measured by Dkhil et alM- 
on a single crystal of PMN, shown as open squares, are 
provided for comparison. 

Several features are interesting to note. First, at low 
temperature the system exhibits an invar-like effect in 
which the cubic lattice parameter changes by less than 
0.001 A; indeed, the data below 320 K are consistent with 
a thermal expansion of zero. At higher temperatures, 
however, the average thermal expansion is 2.5 x 10^^ 1 /K. 
That these disparate regions of null and high rates of 
thermal expansion bracket the revised value for the Burns 
temperature suggests that a direct connection exists be- 
tween the onset of static, short-range polar correlations 
and the structural properties of PMN. This behavior 



seems to be consistent with that of PZN, which also ex- 
hibits an increase in the thermal expansion at temper- 
atures above that where the diffuse scattering first ap- 
pears)^ 

There is ample evidence of similar behavior reported 
by other groups in samples of PMN and PMN-xPT. A 
low-temperature invar-like effect was observed in sin- 
gle crystal PMN-10%PT, where a transition to a high 
rate of thermal expansion was found at 400 K; the ther- 
mal expansion for this sample at high temperature is 
1 X 10~^ 1 /K, which is very close to that measured herei^ 
The x-ray study of ceramic samples of PMN by King et 
al. also show a transition between low and high rates of 
thermal expansion, but larger values for bothf§2 X-ray 
and neutron work conducted by Bonneau et al. on PMN 
powders yielded onset temperatures and values for the 
thermal expansion consistent with our single crystal mea- 
surements.^* Finally, an invar effect was also observed in 
the laser dilatometry study by Arndt and Schmidt on 
a ceramic sample of PMN, which covered a range from 
300 K to 800K.63 

Even though there is a general trend towards a larger 
thermal expansion for temperatures above that where the 
diffuse scattering is onset, there is some sample depen- 
dence and some differences between powders and single 
crystals. As noted by Ye et al.^ powder measurements 
of PMN yield a different slope for the thermal expan- 
sion measurements than do those of Dkhil et al^ Also, 
studies using neutron strain scanning techniques found 
different thermal expansion coefficients as a function of 
depth in single crystal sampleS] ^°'^^ Therefore, part of 
the discrepancy observed between different samples may 
be associated with surface effects. In this regard, we also 
note the presence of a change in the slope of the strain 
curve between 400 K and 500 K; since this change is not 
observed in other PMN studies, we believe this feature 
to be sample-dependent and thus extrinsic. 

The phase-shifted model of polar nanoregions proposed 
by Hirota et al. provides a plausible starting point from 
which to understand the anomalous invar-like behavior 
in PMN below the Burns temperature T^. This model 
is based on the observation that the Pb, Mg, Nb, and O 
ionic displacements obtained from the diffuse scattering 
measurements of Vakhrushev et al^ on PMN can be de- 
composed into a center-of-mass conserving component, 
consistent with the condensation of a soft, transverse op- 
tic mode, and a uniform scalar component, correspond- 
ing to an acoustic phase shift. In so doing Hirota et 
al. were able to reconcile the discrepancies between the 
structure factors of the diffuse scattering and the soft 
TO modes and, in particular, to explain the weakness 
of the diffuse scattering intensity observed in the vicin- 
ity of the (200) Bragg peak^B The idea that the PNR 
are uniformly displaced with respect to the underlying 
cubic lattice in a direction parallel to the PNR polariza- 
tion has already been used by Gehring et al. to explain 
the persistence of the strong diffuse scattering in a single 
crystal of PZN-8%PT against large electric fields applied 
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along [001].^'* In that study it was shown that the diffuse 
scattering near (003), which measures the projections of 
the ionic displacements along [001], decreases under an 
external [001] field as expected. By contrast, the diffuse 
scattering near (300), which measures the projections of 
the ionic displacements along [100], i. e. perpendicular 
to the field direction, remains unaffected even for field 
strengths up to lOkV/cm. Thus a uniform polar state is 
not achieved. This surprising behavior can be understood 
if one assumes that the electric field applied along [001] 
removes or reduces the PNR shifts along [001] while pre- 
serving those along [100]. The diffuse scattering should 
then decrease anisotropically with field, as is observed. 
In the context of the invar-like behavior of PMN below 
Td, we speculate that such uniform shifts of the PNR 
could effectively stabilize the lattice against subsequent 
thermal expansion at lower temperature. Such a scenario 
is consistent with the fact that both PMN-10%PT and 
PMN-20%PT retain cubic structures down to low tem- 
perature when examined with high (/-resolution neutron 
diffraction methodsj ^^'^^ 



VII. DISCUSSION AND CONCLUSIONS 

We have presented a comprehensive neutron study 
of the diffuse scattering in the relaxor PMN. We have 
greatly extended the timescales of previous neutron mea- 
surements by taking data on three different spectrome- 
ters, BT9, SPINS, and the HFBS, which provide elas- 
tic energy resolutions of 500 /zeV, 120 /zeV, and 0.4 /zeV 
HWHM, respectively. While the backscattering data rep- 
resent a 300-fold improvement in energy resolution over 
those obtained by Hiraka et al^ on SPINS, both yield 
an onset temperature of 420±20K for the diffuse scatter- 
ing. This indicates that the PNR in PMN are static on 
timescales of at least 2 ns below 420 K, but are apparently 
dynamic at higher temperatures. Thus the true Burns 
temperature in PMN, which was originally interpreted 
by Burns and Dacol as the condensation temperature 
of static regions of local, randomly-oriented, nanometer- 
scale polarization,^ is 420±20K, not 620 K. Independent 
evidence of the existence of purely dynamic, short-range, 
polar correlations has been reported in PMN-55%PT, a 
composition with no strong diffuse scattering and thus no 
static PNR, by Ko et al. who observed phenomena that 
are typically related to the relaxation of dynamic PNR.^^ 
These include significant Brillouin quasielastic scattering, 
a softening of the longitudinal acoustic phonon mode, 
and a deviation of the dielectric constant from Curie- 
Weiss behavior over an 80 K temperature interval above 

Previous measurements of the temperature depen- 
dence of the neutron diffuse scattering have been ex- 
tended in this study to 900 K, well above Td- In so doing 
we have unambiguously established the existence of two 
types of diffuse scattering based on the observation of 
two markedly different temperature dependencies, one of 



which vanishes at Td and one of which does not. We as- 
sociate the strong, temperature-dependent, diffuse scat- 
tering with the formation of static, short-range, polar 
correlations (PNR) because of its well documented re- 
sponse to an external electric field, and because it first 
appears at the same temperature at which the soft (po- 
lar) TO mode reaches a minimum in energy. We asso- 
ciate the weak, temperature-independent, diffuse scatter- 
ing, which shows no obvious change across either Tc or 
Td, with chemical short-range order because it persists to 
extremely high temperature. We further confirm the ob- 
servations of Hiraka et al.^^ who first characterized the 
distinctly different reciprocal space geometries of both 
types of diffuse scattering, and we show that the bow-tie 
shape of the weak diffuse scattering persists up to 900 K 
»Td. 

Key effects of the strong, temperature-dependent, dif- 
fuse scattering, and thus of the PNR, on the low-energy 
lattice dynamics of PMN are also highlighted in this 
study. The neutron inelastic measurements on PMN- 
60%PT by Stock et al. prove conclusively that PNR can- 
not be the origin of the anomalous broadening of long- 
wavelength TO modes observed in PMN, PZN, and other 
compounds, also known as the waterfall effect, because 
PMN-60%PT exhibits the same effect but no strong dif- 
fuse scattering (no PNR)j2£ By contrast, many stud- 
ies have shown that PNR do broaden the TA modes 
in PMN i4ii5|29 and in PZN-xPT^ but that such ef- 
fects are absent in compositions with no PNR such as 
PMN-60%PT..^ Our cold neutron data show that these 
effects are q-dependent. Whereas long- wavelength TA 
modes with reduced wave vectors q <^ 0.2 remain 
well-defined and exhibit a nearly constant energy width 
(lifetime) from 100 K to 900 K, shorter wavelength TA 
modes with reduced wave vectors q ~ 0.2 A"-'^ broaden 
substantially, with the maximum broadening occurring 
at Td = 420 K. This result motivates a very simple phys- 
ical picture in which only those acoustic phonons having 
wavelengths comparable to the size of the PNR arc signif- 
icantly scattered by the PNR; acoustic modes with wave- 
lengths much larger than the PNR are largely unaffected 
because they simply average over the PNR. Models de- 
scribing this effect have been discussed elsewhere^^S, In 
particular, very recent work by Xu et al. has revealed 
the presence of a phase instability in PZN-4.5%PT that 
is directly induced by such a PNR-TA phonon interac- 
tion. This is shown to produce a pronounced softening 
and broadening of TA modes in those zones where the 
diffuse scattering is strong, and provides a natural expla- 
nation of the enormous piezoelectric coupling in relaxor 
materialsi^ 

In addition, we have performed neutron measurements 
of the thermal expansion with extremely high q- and Tilo- 
resolution over a broad temperature range extending well 
below Tc and far above Td ~ 420 K. In agreement with 
many other studies, our single crystal samples of PMN 
exhibit little or no thermal expansion below T^, behav- 
ior that is reminiscent of the invar effect, but an unusu- 
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ally large thermal expansion coefficient of 2.5 x 10~^ 1/K 
above T^j, where the strong diffuse scattering is absent. 
The crossover between null and large coefficients of ther- 
mal expansion coincides closely with Td, which suggests 
that the appearance of static PNR strongly affects the 
thermal expansion in PMN and thus provides a struc- 
tural signature of the formation of short-range, polar 
correlations in zero field. The model of uniformly dis- 
placed, or phase-shifted, PNR proposed by Hirota et ai, 
which successfully predicts the anisotropic response of 
the strong diffuse scattering to an electric field, offers a 
simplistic, yet plausible, framework in which to under- 
stand this anomalous behavior. 

Finally, it is satisfying to note that the revised value 
of Td = 420 ± 20 K is consistent with the dielectric sus- 
ceptibility data of Viehland et al., from which a Curie- 
Weiss temperature 9 = 398 K was derivedi^^ Such good 
agreement between Td and O solidifies our identification 
of the strong diffuse scattering with the condensation of 
the soft TO mode, which reaches a minimum frequency 
at Td = 420 K. However this begs the question of how 
one should interpret the original Burns temperature of 
^ 620 K. At present there are two broadly divergent opin- 
ions on this issue, one of which considers ~ 620 K to be a 
meaningful temperature scale in PMN, and one of which 
does not. As it turns out, this debate is closely tied to 
another on how many temperature scales are needed to 
describe the physics of relaxors. We offer no final res- 
olution to this discussion. Instead, we close our paper 
with a brief summary of the primary studies supporting 
these contrasting points of view, which is by no means 
comprehensive, so that the readers may draw their own 
conclusions. 

A number of experimental studies report evidence of 
either structural or dynamic changes in PMN in the tem- 
perature range 600 K - 650 K, starting with the opti- 
cal index of refraction measurements of Burns and Da- 
col.^ Siny and Smirnova were the first to observe strong, 
first-order Raman scattering in PMN at high tempera- 
tures, which, being forbidden in centrosymmetric crys- 
tals, implied the presence of some type of lattice dis- 
tortion.^^ In 2002 Vakhrushev and Okuneva calculated 
the probability density function for the Pb ion in PMN 
within the framework of the spherical layer model us- 
ing x-ray powder diffraction datai^I It was shown that 
this probability density evolves from a single Gaussian 
function centered on the perovskite A-site to a double 
Gaussian form between 635 K and 573 K, and that the 
positions of maximum density for the lead ion follow 
a power law {Td - Tf '^^ with Td = 635 K. This pic- 
ture was developed further by Prosandeev et al. within 
a model that ascribed the changes in the lead probabil- 
ity density function and subtle variations in the thermal 
expansion near 620 K to a crossover from soft-mode to 
order-disorder dynamics;^ In 2004 Egami et al. first pro- 
posed that at very high temperatures the vibrations of 
the oxygen octahedra are sufficiently faster than those 
of the Pb ions that, on average, a high-symmetry, local 



environment is obtained, whereas at temperatures below 
~ 680 K the Pb and O ions become dynamically corre- 
latedi^^ This led to the subsequent reinterpretation of the 
Burns temperature (~ 600 K) as being the point below 
which dynamic PNR first form, which was based on the 
dynamic PDF measurements of Dmowski et al. obtained 
with an elastic energy resolution of 4 meV at six temper- 
atures (680 K, 590 K, 450 K, 300 K, 230 K, and 35 K).^ 
Interestingly, after integrating their data from -5 meV to 
-f 5 meV, Dmowski et al. saw evidence of a third tem- 
perature scale in PMN of order 300 K, and thus distinct 
from both T^ ~ 210 K and Td ~ 600 K, which they as- 
sociated with the formation of static PNR. Similar ideas 
have very recently been put forth by Dkhil et al. based 
on extremely interesting acoustic emission and thermal 
expansion measurements, and also in a general summary 
written by Toulouse . ^^'^^ 

A different approach was taken by Stock et al. who 
proposed a unified description of the lead-based relax- 
ors on the basis of striking similarities between PMN 
and the Zn-analogue PZN, both of which display strong, 
temperature-dependent, diffuse scattering; identical soft 
mode dynamics; yet no long-range, structural distortion 
at low temperature in zero field. ^^'^^ Arguing by analogy 
with magnetic systems, Stock et al. considered a three- 
dimensional Heisenberg model with cubic anisotropy in 
the presence of random fields in which the Heisenberg 
spin corresponds to the local ferroelectric polarization, 
the cubic anisotropy represents the preferential orien- 
tation of the polarization along the axes, and 
the isotropic random magnetic field corresponds to the 
randomly-oriented local electric fields that originate from 
the varying charge of the B-site cation. Following a 
suggestion by Aharony, Stock et al. considered the case 
where the Heisenberg term dominates over the random 
field term, and the cubic anisotropy term is the weak- 
est of the three. In this picture there would be just two 
distinct, static, temperature scales. For T > Td, the 
cubic anisotropy is irrelevant and therefore the system 
should behave like a Heisenberg system in a random field. 
In this case the excitation spectrum is characterized by 
Goldstone modes and therefore no long-range order is 
expected in the presence of random fieldsi^ Instead the 
system forms polar nanoregions in a paraelectric back- 
ground. The second temperature scale Tc appears at low 
temperatures where the cubic anisotropy term becomes 
important, and in this limit the system should resem- 
ble an Ising system in the presence of a random field. 
This model thus explains the local ferroelectric distor- 
tion characterized by the recovery of the soft-optic mode 
and, although a 3D Ising system in equilibrium should 
display long-range order in the presence of small random 
fields, as is observed in magnetic systems, nonequilibrium 
effects with long time scales become dominant. The pres- 
ence of such nonequilibrium effects may explain the lack 
of long-range ferroelectric order in PMN and PZN as well 
as the history dependence of physical properties such as 
the linear birefringence. The phase-shifted nature of the 
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polar nanoregions may also create another energy barrier 
against the ordered phase at T^. Stock et al. are also able 
to explain the physics of compounds beyond the MPB, 
such as PMN-60%PT, within this model, for which only 
one temperature scale (Tc) exists. As the Ti concentra- 
tion Ti increases, the relaxor phase diagram crosses over 
to a ferroelectric regime, and this can be understood as an 
increase in the strength of the cubic anisotropy term that 
is simultaneously accompanied by a reduction of the ran- 
dom fields as must occur in the limit of pure PbTiOs.- - 
It should be mentioned that other studies have invoked 
random fields to explain the static and dynamic proper- 
ties of relaxors:5M2,i22, 

Independent of the validity of either of the two pic- 
tures summarized above, the seminal finding of our study 
of remains that the strong diffuse scattering in PMN 
first appears at a temperature that depends sensitively 
on the energy resolution of the measurement. This fact 
inevitably raises interesting questions about the signif- 
icance of the previous value of the Burns temperature 
(~ 620 K). If the strong diffuse scattering in PMN results 
from the soft TO mode, then other techniques based on 
x-ray diffraction, thermal neutron scattering, or neutron 
PDF, which provide comparatively much coarser energy 
resolution, will be unable to discriminate between low- 
energy, dynamic, short-range polar correlations and truly 
static {fiLo — 0) PNR because any low-energy, polar cor- 
relations will be integrated into the elastic channel by 
the broad energy resolution. Thus as the TO mode soft- 
ens, the associated low-energy, polar, spectral weight will 
fall within the energy resolution at temperatures above 
Td = 420 K, and the net result will be an artificially 
higher value of Td that increases with the size of the 
energy resolution. This effect should be especially pro- 
nounced for phonon modes that are broad in energy (i. e. 
that have short lifetimes), and this is clearly the case for 
the soft TO mode in PMN, which exhibits a nearly over- 
damped lineshape and a minimum frequency precisely at 
Td = 420 K. While it is true that the structure factor 
of the strong diffuse scattering is inconsistent with that 
of the soft TO mode, the phase-shifted model of Hirota 
et al. provides a way to reconcile this discrepancy. In 
this respect, we simply suggest that the previous value of 
Td ~ 620 K might not represent a physically meaningful 



temperature scale in PMN. Future studies examining the 
values of Td in other relaxor compounds with improved 
energy resolution are thus of interest. 

Our reassessment of the Burns temperature Td imme- 
diately clarifies an intimate relationship between the for- 
mation of static, short-range polar correlations and the 
consequent effects on both the low-energy lattice dynam- 
ics and structure of PMN. Cold neutron triple-axis and 
backscattering spectroscopic methods conclusively show 
the existence of static, short-range polar correlations, 
only below Td = 420 ± 20, K on timescales of at least 
2 ns. Thermal neutron measurements of the lattice dy- 
namics reflect the presence of these static PNR through 
the presence of a distinct minimum in both the soft TO 
and TA mode energies, both of which occur at 420 K. The 
effect of PNR on the lattice dynamics is evident only for 
TA modes having wave vectors of order 0.2 A~^, a fact 
that could be exploited to determine the size of the PNR. 
At the same time an enormous change in the coefficient 
of thermal expansion is seen near Td, below which the 
crystal lattice exhibits invar-like behavior. 
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